Anisotropic scaling and generalized conformal invariance at Lifshitz points 
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The behaviour of the 3D axial next-nearest neighbour Ising (ANNNI) model at the uniaxial Lifshitz 
point is studied using Monte Carlo techniques. A new variant of the Wolff cluster algorithm permits 
the analysis of systems far larger than in previous studies. The Lifshitz point critical exponents are 
a = 0.18(2), P = 0.238(5) and 7 = 1.36(3). Data for the spin-spin correlation function are shown to 
be consistent with the explicit scaling function derived from the assumption of local scale invariance, 
which is a generalization of conformal invariance to the anisotropic scaling at the Lifshitz point. 
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The modern understanding of critical jphenomena is 
governed by the notion of scale invarianceE For isotropic 
critical systems the extension from global, spatially ho- 
mogeneous scaling to space-dependent rescaling factors 
leads to the requirement of conformal invariance of the 
correlators. This approach has proven to be very fruitful 
when examining isotropic eojii^ibrium critical systems, 
especially in two dimensions .EJu 

On the other hand, little is known for systems with 
strongly anisotropic critical points, where the value of the 
anisotropy exponent 9 differs from unity. In these cases, 
the two-point function C{fj_, rn ) satishes the scaling form 



C(fL,r||) = r2^C(6fL,6^||) =rj2^0(r||r] 



(1) 



where and rj_ ~ |fj_| are the distances parallel and 
perpendicular with respect to a chosen axis, a; is a scaling 
dimension, 9 is the anisotropy exponent and f2(u) is a 
scaling function. Scale invariance alone is not enough to 
determine the form of the function ri(w). 

Recently, a generalization of conformal invariance in- 
volving local space-time-dependent scale transforniations 
for anisotropy exponents 9^1 has been proposedJj This 
approach attempts to generalize the scaling (|l|), usually 
considered with h constant, to space-dependent rescaling 
b = b{f_\_, ry), thereby assuming that the two-point func- 
tions still transform in a simple way. These transforma- 
tions are constructed, starting from the conformal trans- 
formations ry (ary +/3)/(7ry +6) with aS — Pj = 1, in 
such a way that the transformations in the 'spatial' co- 
ordinates f± are consistent with the anisotropic scaling 
(^ . Systems which are invariant under these transforma- 
tions and whose correlators, generalizing (^), transform 
covariantly under them, are said to satisfy local scale in- 
variance (LSI). It turns out that if 9 — 2/N, where N 
a positive integer, f2(w) satisfies the differential equatio: 
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where ( — 2x/9 aad ai is a constant. Eq. (g) can 
be explictly solvedu in terms of hypergeometric func- 
tions 2PN-1 (conformal invariance is reproducedo for 



N = 2 and = 1 gives the non-relativistic Schrodinger 
invarianceO) . Evidently, the above hypothesis of LSI in 
systems satisfying (|^) is a strong one and relies on cer- 
tain assumptions about the structure of the underlying 
field theory. In this letter, we shall test the idea of LSI 
by explicitly checking the resulting expressions for il{v) 
in a non-trivial spin system which satisfies the strongly 
anisotropic scaling (0). 

While dynamical scaling (^ occurs in critical dynamics 
(then 9 is referred to as dynamical exponent) or in true 
non-equilibrium phase transitions such as directed perco- 
lation, well-known examples of stongly apisotropic equi- 
librium criticality are the Lifshitz pointM encountered in 
systems with competing interactions. At a Lifshitz point, 
a disordered, an uniformly ordered and|-|a periodically 
ordered phase become indistinguishable.Q The simplest 
model for these is the ANNNI (axial next-nearest neigh- 
bour Ising) modelEl which describes faithfully, among 
others, magnetic systems, alloys or uniaxially modu- 
lated ferroelectrics.BT2l Recently, a large varieto. of new 
physical systeir^. (ferroelectric liquid, crystals,El uniax- 
ial ferrolectrics,L3 block copolymersEj or even quantum 
systemstil) were shown to possess a Lifshitz point which 
has stimulated renewed interest p. its properties. Fur- 
thermore, new field theory studiesEj~t3 have lead to more 
refined estimates (in the framework of an e-expansion) 
of the critical exponents of the general m-fold Lif- 
shitz pointp-|iii-|rf dimensions with a rt-component order 
parameter .EjiiZI 

Here, we study the 3D ANNNI model, defined on a 
cubic lattice with periodic boundary conditions. The 
Hamiltonian is 

W = Sxyz {s(x+l)yz + Sx(y+l)z + Sxy(z+1)) 

xyz 

(3) 

xyz 

with Sxyz = ±lj whereas J > and k > are coupling 
constants. In z-direction competition between ferro- 
magnetic nearest neighbour and antiferromagnetic next- 
nearest neighbour couplings takes place, leading to a 
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rich phase diagram with a multitude of different phases, 
both commensurate and incommensurate to the under- 
lying cubic lattice.Q The anisotropy exponent 9 = v^/v^, 
where i^y und are the exponents of the correlation 
lengths parallel and perpendicular to the z-axis. At the 
uniaxial Lifshitz point, a recent careful field-theoretical 
calculationllB shows that 6 = \ — ae^ -\-0{e^) in a second- 
order e-expansion (with e = 4.5 — d) where a ~ 0.0054 
for the 3D ANNNI model. 

Our main purpose will be the numerical computation 
and thorough analysis of the critical spin-spin correla- 
tion function at the uniaxial Lifshitz point of the ANNNI 
model through a large-scale Monte Carlo simulation. The 
agreement between our numerical results and the ana- 
lytic expression for Q,{y) derived from (H) presents evi- 
dence that local scale invariance, as formulated in Ref. 

is realized as a new symmetry in strongly anisotropic 
equilibrium critical systems. 

Such a study does require reliable and precise estimates 
of the critical exponents. However, published estimates 
of critical exponents obtained with different techniques 
spread considerably, see Table I. We therefore undertook 
large-scale Monte Carlo simulations to estimate, the ex- 
ponents reliably. Previous Monte Carlo studieaij consid- 
ered only small systems of (mostly) cubic shape. Here, 
we present calculations for large systems of anisotropic 
shape with L x L x N spins, with 20 < L < 240 and 
10 < < 100, taking into account the special finite-size 
effects-coming from the anisotropic scaling at the Lifshitz 
point nil This is the first study of the ANNNI model where 
the exponents a, /3, and 7 are computed independently. 

As usual, the problems coming from critical slowing- 
down encountered when using local Monte Carlo dynam- 
ics, are alleviated by using _npn- local methods, such as 
the Wolff cluster algorithm^ For the Ising model with 
only a nearest neighbour coupling J, this algorithm may 
be described as follows: one chooses randomly a lat- 
tice site, the seed, and then builds up iteratively a clus- 
ter by including a lattice site j (with spin Sj), neigh- 
bour to a cluster site i (with spin s^), with probability 
p = 1(1 + sign (siSj)) (1 - exp [-2J/{kBT)]). One ends 
up with a cluster of spins having all the same sign which 
is then flipped as a whole. This kind of same-sign clus- 
ters are obviously not adapted to our problem because of 
the competing interactions along the z-direction, see (|^). 

We therefore propose the following modified cluster 
algorithm. Starting with a randomly chosen seed, one 
again builds up iteratively a cluster. Consider a newly 
added cluster lattice site i with spin Si. A lattice site 
j with spin Sj nearest neighbour to i is included with 
probability Pn = p, whereas an axial next-nearest neigh- 
bour site k with spin Sk is included with probability 
Pa = i(l- sign (s,Sfe))(l- exp [-2jK/(fcsr)]). Thus, 
the final cluster, which will be flipped as a whole, con- 
tains spins of both signs. Ergodicity and detailed balance 
are proven as usual. This algorithm works extremely well 
in the ferromagnetic phase and in the vicinity of the Lif- 
shitz point. Generalization to other systems with com- 



peting interactions is straightforward. For the computa- 
tion of the spin-spin correlation function we adapt in a 
similar way a recently proposed very efficient algorithm 
using Wolff clusters .lJ This algorithm yields the infinite- 
system correlation functions at temperatures above 
and largely reduces finite-size effects at the critical tem- 
perature as compared to a more traditional approach. 
These algorithms will be discussed in detail elsewhere.a 
We now outline the determination of the critical expo- 
nents. The results are listed in Table I. As an example, 
Figure 1 shows the effective exponent 



where m denotes the magnetization and t = 1 — T/Tc. 
In the limit the effective exponent yields the crit- 
ical exponents (3, provided finite-size effects can be ne- 
glected. The two sets of data in Figure 1 correspond to 
two different paths in the temperature-interaction space, 
both ending at the point (k = 0.270, Tc = 3.7475), set- 
ting J/fcs = 1. For set (a) k was fixed at 0.270, whereas 
for set (b) n = 0.270-^1.6 (l/T-l/rj. The corrections to 
scaling for set (b) are small compared to set (a), resulting 
in a plateau for t < 0.06, thus making a very precise esti- 
mation of (3 possible. Of course, finite-size effects have to 
be monitored carefully. As usual, we adjust the systeHi 
size in order to circumvent finite-size dependencesEj'E^I 
For the determination of the susceptibility and specific 
heat critical exponents 7 and a, see Table I, data ob- 
tained at temperatures both below and above Tc were 
analysed. Our error bars take into account the sample 
averaging as well as the uncertainty in the location of the 
Lifshitz point. Based on our data, we locate the Lifshitz 
point at K = 0.270 ± 0.004, = 3.7475 ± 0.0005, thus 
confirmina-an estimation from a high temperature series 
expansion.Ea 

The agreement of the independently estimated expo- 
nents a, (3 and 7 with the scaling relation a-\-2(3-\-^ — 2, 
up to « 0.8%, illustrates the reliability of our data. 

We are now ready to discuss the scaling of the spin- 
spin correlation function C(7tL,7'j|) = (^sp^^r\\ ^00) ^^"^ 
its scaling function n{v) as defined in (|l]). In (d_L -I- 1) di- 
mensions, one has C = 2(d_L -t- 0)/e{2 + 7//?). For the iD 
ANNNI model, C = 1-30 ± 0.05, where the error follows 
from the errors in the values of the critical exponents. In 
Figure 2 we show selected data for the function 

^{u)^u-^Q.{l/u) (5) 

with u — /ry , as computed by Monte Carlo simul^ 
tions of a system with 200 x 200 x 100 spins (assumingEZi 
9 — 1/2). This permits a nice visual test of the data 
collapse and establishes scaling. 

The small deviations (of order « 2%) from the value 
6* = 1/2 obtained in recent field-theoretical calculationst3 
are not yet distinguishableE3 from the purely numerical 
errors in our data and the exponents derived from them. 
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Therefore, for our purposes, namely the test of LSI, it 
is enough to set 6 = 2/N = 1/2, leading to = 4 
in the differential equation (^). In addition, the scaling 
form implies the boundary condition Q{v) ^ v~'' for 
i; ^ oo. For = 4, there are two indepcndeiit-Bolutions 
of cq. (H) satisfying this boundary conditionB'Q and the 
scaling function becomes 



n{v) = 5oi^o(«/(4ai)'/') + 6i«i^i(z;/(4ai)i/4) (6) 



whereu 



r(C/4)^n!r(n/2 + 3/4) 



Fiix) 



r(3/2) 



I r((n + l + C)/4)s(n) 

r(i(C + i))^or(n/4 + i)r(i(n + 3)) 



(7) 



and s{n) = -^(cos^ + sin if) cos if. Here ai is 

the constant occuring in and 6o,i are free parame- 
ters. Since bo and ai are merely scale factors, the func- 
tional form of the scaling functions n{v) and $ (u) — 
n{l/u) depends on the single universal parameter 



-C 



1/4 



bi/bo. 



To see this, consider the moments M{n) :— 
Jl^ duu"^{u). For 6* = 1/2 and taking into account (|^), 
it followf£3 that in the scaling region the moment ratios 



k Ik 

i=i / i=i 



(8) 



with k > 2 and J2i = J2j ^^'^ independent of 5o 
and ai and only depend on the functional form of 
as parametrized by p. Our Monte Carlo data for the spin- 
spin correlator will be consistent with LSI if the values 
of p determined from several different ratios J coincide. 

As we are not able to compute numerically the func- 
tion for values of u below uq « 0.22 a direct 
analysis along the lines just sketched is not possible. 
Instead we have to consider the moments M{n) := 
/„'^duu"$(ua}/^) = a^^''+^^'^ Qdww^^iw) with 

wq = uoa}^^. The moment ratios J({mi} ; {nj}) (de- 
fined in complete analogy with (||)) then depend on the 

1 /4 

scale factor through wq. The parameters ai and 
p are determined from the following scheme. Choosing 
a suitable starting value for ai we compute an approx- 
imative value for p by comparing the values of the mo- 
ment ratios J (obtained from the full data set for <I>(w) 
as shown in the inset of Figure 2) with the p-dependent 
expressions coming from integrating the M(n) using the 
analytic form (^,|^. An improved value for ai is then de- 
rived by comparing the values M{m)/ M{n) for arbitrary 
m and n obtained (i) from our numerical data and (ii) 
from the analytical expressions after inserting the value 
of p. The final values of ai are obtained by averaging 
over five different pairs (m, n). 



The values of p and ai determined from several distinct 
moment ratios are collected in Table II. We obtain the 
mean values p — —0.11(1) and ai = 33.2(8). The consis- 
tency of the different determinations of the two parame- 
ters provides clear evidence in favour of the applicability 
of the hypothesis of local scale invariance to the Lifshitz 
point of the ANNNI model. 

Finally, foo ~ 0.41 is obtained by adjusting the scale of 
$. Inserting these values into the analytical expression 
yields for <& the full curve shown in the inset of Figure 2. 
The agreement between our MC data and the theoretical 
result is remarkable. 

Local scale invariance had been confirmed before at 
the Lifshitz point of the ANNNS model. In that ex- 
actly solvable model, the Ising model spins in (^ are re- 
placed by spherical model spias Sxyz € IR together with 
the usual spherical constraint!] At the Lifshitz point, opa 
has 6 = 1/2 and the exactly known spin-spia correlatoicJ 
agrees with the scaling form (^ for hi = O.Q Our finding 
that LSI also appears to hold for the ANNNI model sug- 
gests that the domain of validity of LSI should extend 
beyond the context of free field theory which underlies 
the ANNNS model. It appears plausible that LSI will 
also hold true for the Lifshitz points of the ANNNXY, 
ANNNH,. . . modelsQ which are intermediate between the 
ANNNI and the ANNNS model. Since the number of 
dimensions d±^ merely enters as a parameter, local scale 
invariance pGOuld also be tested along the lines of an e- 
expansion.EJ Finally, the consistency of other correlators 
(e.g. energy-energy) with LSI should be tested. We plan 
to come back to this elsewhere. It would be tempting 
to see whethet, the powerful techniques of 2D confor- 
mal invarianccD might be extended to the situation of 
anisotropic scaling realized at Lifshitz points. This would 
lead to numerous physical applicationsH llj and is under 
investigation. 

In conclusion, the precise localisation of the 3-D 
ANNNI model Lifshitz point and improved estimates of 
its critical exponents allowed for the first time to deter- 
mine reliably the scaling of the spin-spin correlator. Its 
functional form was found to agree with the prediction of 
local scale invariance. The confirmation of the applica- 
bility of local scale invariance to this situation suggests a 
new symmetry principle for the description of equilibrium 
systems with anisotropic scaling, especially for systems 
with competing interactions at their Lifshitz points. 
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TABLE L Critical exponents at the Lifshitz point of the 
3D ANNNI model, as obtained from Monte Carlo simulations 
(MC), high-temperature series expansion (HT) and renormal- 
ized field theory (FT). The numbers in brackets give the es- 
timated error in the last digit (s). 
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1.315 
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0.18(2) 


0.238(5) 


1.36(3) 


1.34(5) 


0.175(8) 



TABLE II. Values of the parameters p and ai computed 
from different moment ratios J{{mi} ; {nj}), see text. 
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FIG. 1. Effective exponent /3off versus t for two different 
trajectories in the {T,n) space, see text. Error bars include 
the uncertainty in Td Tc{n = 0.270) = 3.7475 ± 0.0005. 
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Scaling funtion $(u) (see text) 
(yrT/rii)^/^ for K = 0.270 and T = 3.7475. 
Selected Monte Carlo data for a system of 200 x 200 x 100 
spins are shown.EZI The different symbols correspond to the 
values of r±. Inset: comparison of the full data set of 1.7- 10* 
points for the scaling function $(«) (gray points) with the 
analytical prediction eqs. (^J^ following from the assumption 
of LSI, with p = -0.11, ai = 33.2 and bo = 0.41 (full curve). 
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